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SECTION  1 


GENERAL  DESCRIPTION 


1 . 1  Purpose  of  the  Program  Maintenance  Manual 

The  objective  of  this  Program  Maintenance  Manual 
for  the  Dynamic  Thermodynamic  Sea  Ice  Model,  PIFS-N,  is  to 
provide  the  maintenance  programmer  personnel  with  the  information 
necessary  to  effectively  maintain  the  system. 


1 . 2  Background 

Sea  ice  forecasting  programs  have  been  clearly 
connected  with  fleet  polar  operations  since  the  early  1950's. 
When  organizing  and  conducting  the  sea  lift  operations  reguired 
for  the  establishment  and  resupply  of  arctic  bases  such  as  Thule, 
Greenland  and  the  Distant  Early  Warning  Line  across  Alaska  and 
Canada,  ice  reconnaissance  and  forecasting  services  were 
requested . 


Since  the  advent  of  the  underice  submarine  operations  in 
1957,  ice  intelligence  prediction  services  have  been  requested  by 
suDmarine  commands.  Their  operational  forecast  requirements 
deviated  markedly  from  surface  ship  or  icebreaker  requirements. 
The  submarine  operator  is  primarily  interested  in  knowing  the 
distribution,  in  frequency  and  size,  and  depth  of  ice  pressure 
ridges  which  may  constitute  a  hazard  to  underice  navigation 
particularly  in  shoal  water. 

The  surface  ship  is  primarily  interested  in  knowledge  of 

I  ice  concentration,  the  distribution  of  the  various  stages  of  ice 

development  and  floe  sizes. 

|  Escalations  of  operational  activities 
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in  ice-covered  waters  over  the  past  seven  years  with  attendant 
navigation  problems  have  focused  Navy  and  national  attention  on 
the  need  for  a  more  reliable  sea  ice  forecasting  program. 

The  recent  establishment  of  the  Navy  NOAA  Joint 
Ice  Center  is  a  strong  indication  that  Navy  and  other 
governmental  departments  are  moving  forward  to  meet  existing  and 
full  operational  requirements. 

In  recent  years,  important  contributions  have  been 
made  to  understand  the  dynamics  of  sea  ice  through  the  Arctic  Ice 
Dynamics  Joint  Experiment  (AIDJEX) .  Several  mathematical  models 
have  been  developed.  A  continuing  effort  is  needed  to  evaluate 
the  applicability  of  these  models  for  sea  ice  forecasting. 

The  dynamic-thermodynamic  Sea  Ice  Model  was 
developed  by  the  W.  D.  Hibler  III  (Hibler,  1979).  The  Naval  Air 
Systems  Command  (NAIR-270G)  tasked  the  Naval  Ocean  Research  and 
Development  Activity's  Polar  Oceanography  Branch  to  implement  the 
model  utilizing  the  Fleet  Numerical  Oceanography  Center  (FNOC) 
environmental  data  base  and  Cyber  203. 

The  model  uses  atmospheric  forecast  and  analysis 
data  available  at  FNOC  in  the  opeBltional  data  base.  The  model 
outputs  forecasts  of  ice  drift,  concentration,  thickness, 
convergence/divergence  plus  ice  and  open  water  growth. 
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There  have  been  a  number  of  sea  ice 
models  appearing  in  the  literature  in  the  recent  years. 
A  list  of  applicable  references  to  this  model  is  presented 
below: 

i)  Hibler,W.D.,  1980:  Modeling  a 

Variable  Thickness  Sea  Ice  Cover.  Mon ,  Wea.  Rev. ,  108, 

1943-1973. 

ii)  Hibler,  W.D.  1979:  A  dynamic  sea 
Ice  Model.  J.  Phys.  Oceanoqr.,  9,  815-846. 

iii)  Pritchard,  R.S.,  1978:  The 
effect  of  strength  on  simulation  of  sea  ice  dynamics. 
Proc.  Fourth  int.  Conf.  Port  and  Ocean  Engineerin  under 
Arctic  conditions,  D.E.  Muggeridge,  Ed.,  Memorial 
University  of  St.  Johns,  Newfoundland,  494-505. 

iv)  Pritchard,  R.S.,  M.D.  Coon  and  M.G. 

McPhee,  1977 :  Simulation  of  sea  ice  dynamics  during 

AIDJEX.  J.  Prep.  Vessel  Tech.,  99J,  491-497. 

v)  Thorndike,  A.S.,  and  R.  Colony, 

1980:  Large-scale  ice  motion  in  the  Beaufort  Sea  during 

AIDJEX,  April  1975  -  April  1976.  Sea  Ice  Processes  and 
Models,  R.S.  Prichard,  Ed.,  University  of  Washington 
Press,  249-260. 

vi)  FNWC  User  Manual 

vi)  Computer  Operational  Manual. 

1.3 

FNOC 
NORDA 

NPOC 
NEDN 


Terms  and  Abbreviations 

Fleet  Numerical  Oceanography  Center 
Naval  Ocean  Research  and  Development 
Activity 

Naval  Polar  Oceanography  Center 
Naval  Environmental  Data  Network 


Section  2 


SYSTEM  DESCRIPTION 


2.1  System  Application 

Fleet  Numerical  Oceanography  Center  is 
a  large  computer  complex,  tasked  with  the  mission  of  pro¬ 
viding  worldwide  environmental  support  to  the  U.S. 
fleets.  To  accomplish  this  mission,  FNOC  collects  meteo- 
rogical  report  data,  analyzes  the  data  and  predicts  changes 

in  environmental  conditions. 

The  Dynamic  Thermodynamic  Sea  Ice  model 

developed  by  Hibler  is  modified  to  use  FNOC  environmental 
data.  The  modified  model  gets  input  data  from  the  FNOC 
data  base.  These  data  consist  of  wind  data,  surface 
pressure,  surface  vapor  pressure,  air  specific  humidity, 
air  temperature,  incoming  and  outgoing  radiation  on  the 
FNOC  63  x  63  Northern  Hemisphere  grid  (A01,A10,  All,  A12, 
A16,  A16,  A20 ,  A21)  ,  interpolates  these  data  to  the  model 
grid  to  get  initial  boundary  conditios. 

The  model  computes  ice  drift,  ice  con¬ 
centration,  thickness,  convergence/divergence  rate  plus 
thick  and  thin  ice  growth/decay  rate. 


2. 2  Security  and  Privacy 

The  sea  ice  model  does  not  currently 
have  access  to  classified  information  nor  does  it  produce 
classified  output.  Hence,  security  should  be  treated  at 
appropriate  levels,  and  the  user  is  responsible  for  pro¬ 
tection  of  any  material  used  by  the  model. 


2.3 


The  dynamic  thermodynamic  sea  ice  model 
has  been  modified  to  run  on  the  Cyber  203  utilizing  the 
FNOC  environmental  data  base,  FNOC  library  and  plotting 
facility.  In  general,  the  modified  model  can  be  divided 
into  three  modules:  -Input  Module,  Processing  Module, 

Output  Module. 

Input  module  creates  a  model  grid  on 
the  FNOC  63  x  63  Polar  Stereographic  grid,  sets  boundary 
and  initial  conditions  for  the  model  and  obtains  and  in¬ 
terpolates  input  data  to  each  model  grid  point. 

The  processing  module  can  be  divided 
into  two  parts  which  process  two  different  mechanisms: 
the  dynamic  mechanism  and  the  thermodynamic  mechanism. 

The  heart  of  the  dynamic  part  is  sub¬ 
routine  RELAX  which  computes  ice  drift  for  the  model. 

The  thermodynamic  mechanism  is  pro¬ 
cessed  by  subroutine  HEAT  which  estimates  ice  thickness, 
concentration  and  thick  and  thin  ice  growth/decay  rate. 

The  output  module  calls  many  FNGC 
routines  to  format  output  data  into  CRANDIC  format,  and 
stores  them  for  plotting  and  printing.  A  functional  dia¬ 
gram  of  the  model  is  shown  in  Figure  1.  Figure  2 
illustrates  more  detail  of  each  functional  division. 
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Figure  1 


The  functional  diagram  of  the 
Dynamic  Thermodynamic  Sea  Ice  Model 


Input  Module 
Create  the  model  grid 
Set  up  boundary  and  initial  conditions 
Get  input  data 


- * _ 

Processing  Module 
Dynamic  Part: 
Calculate  ice  drift 
Thermodynamic  Part: 
Calculate  ice  thickness, 
coverage,  growth/decay  rate 


_ ± _ 

Output  Module 
Store  output  on  CRANDIO 
for  plotting  and  print 
output 


Figure  2 


The  chart  of  interrelationship  of  the 
major  components  of  the  system 


Create  the  model  grid 
Set  up  boundary  and  initial  conditions 
Get  input  data:  winds,  dynamic  height, 
currents,  surface  pressure 
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Estimation  of 

ice 

evolution  of  ice 

strength  from 

the  ice 

thickness  charcater- 

thickness  and 

concen- 

istics 

tration 

Thermodynamic  Model 


Thermodynamic  input: 
air  tempreature , solar 
radiation,  air  pressure 


Output:  ice  thickness, 

drift,  coverage,  ice 
growth  rate,  stored  in 
ZRANDIO.  and  printed 
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2.4  Program  Description 

The  following  list  of  subroutine  names 
and  functions  is  provided  for  easy  reference: 


NAME 

PROGRAM  ICEMDL 


FUNCTION 

Main  driving  program  for  the 
model 


SUBROUTINE  MESH 

BNDRY 

OCEAN 


Input  Module 

Calculates  FNOC  (I,  J)  grid 
points  for  the  model  grid; 
Sets  up  boundary  masks; 
Computes  ocean  surface 
currents; 


INITIAL  Obtains  initial  atmospheric 

data. 


SUBROUTINE  FORM 


XSUM 

PLAST 


UVDDFF 

INTRP 


Procesing  Module 

Sets  up  forces,  drag  co¬ 
efficients,  non-linear 
viscosities  for  use  in  each 
time  step; 

Sums  a  vector; 

Calculates  non-linear  vis¬ 
cosities  for  plastic  flow, 
used  in  FORM 

Converts  U,V  to  direction 
and  magnitude  and  vice  versa 
Performs  interpolation  from 
the  FNOC  grid  to  the  model 
grid; 
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RELAX 


SUBROUTINE 


FELLIP 

FELLDI 

ADJUST 

MEAN 


ADVECT 

DIFFUS 

HEAT 

BUDGET 

AVG 


GROWTH 


DIVERG 

UVPLOT 

HAPLOT 

PRNT 

STATPRT 

GROWPEC 


Solves  linearized  momentum 
balance  with  spatially  varying 
bulk  and  shear  viscosities. 
Uses  sequential  over-re¬ 
laxation  technique; 

Calculates  finite  differences 
for  use  in  RELAX; 

Estimates  thickness  and  con¬ 
centration  at  "open  boundary" 
grid  cells  based  on  adjacent 
values ; 

Performs  advection  of  ice 
thickness  and  concentration; 
Driving  subroutine  for  thermo¬ 
dynamic  calculations; 

Computes  thin  and  thick  ice 
growth/decay  rate; 

Computes  averaged  input 
variables  for  staggered  grid 
system; 

Calculates  changes  in 
concentration  and  mean  ice 
thickness  due  to  growth  and 
redistribution  due  to  ridging. 

Output  Module 

Output  convergence/divergence 
in  CRANDIO  format; 

Outputs  direction  and  magni¬ 
tude  of  ice  drift  in  CRANDIC 
format ; 

Outputs  ice  thickness  and  con¬ 
centration  in  CRANDIO  format 
Provides  hard  copy  print  out 
of  model  variable  arrays; 

Prints  resource  use 
statistics; 

Outputs  ice  growth/decay  rates 
in  CRANDIO  format. 
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2.4.1 


*1 


Program  ICEMDL 

A  brief  description  of  the  theory  in¬ 
volved  in  this  model  is  initially  presented  in  order  to 
help  clarify  and  provide  a  reference  for  the  subroutine 
descriptions  subsequently  presented. 

The  main  driving  program,  ICEMDL,  is 
designed  to  call  the  subroutines  which  set  up  the  initial 
conditions,  perform  the  equation  integration  and  output 
the  specified  forecasts  in  CRANDIC  format  on  the  CDC  CYBER 
203  at  FNOC. 

The  overall  structure  essentially  con¬ 
sists  of  three  main  components: 

i)  input  and  initial  conditions; 

ii)  equation  integration; 

iii)  output  processing. 

The  initial  conditions  are  determined 
from  the  atmospheric  data  available  through  the  FNGC  oper¬ 
ational  data  base.  A  model  grid  is  defined  to  cover  the 
arctic  ocean  basin  and  is  a  subset  of  the  FNGC  hemispheric 
grid  (Figure  1)  .  the  grid  contains  a  square  mesh  with  a 
distance  of  approximately  120  Km  between  grid  points. 
Atmospheric  variables  are  interpolated  from  the  FNOC 
hemispheric  grid  to  the  model  grid  through  the  use  of  a  16 
point  Bessel  interpolation  scheme.  The  following  atmos¬ 
pheric  parameters  are  interpolated  to  the  model  grid; 


i) 

surface 

wind  (u,v) ; 

ii) 

surface 

pressure ; 

iii) 

sur  face 

vapor  pressure; 

iv) 

surface 

air  temperature; 

v) 

incoming 

solar  radiation; 

vi) 

total  heat  flux  at  the  surface 

vii) 

sensible 

plus  evaporative  heat 

flux  at 

the  surface. 

Input  processing  begins  by  determining 
the  boundary  masks  to  be  used  in  the  simulation.  The 
boundary  masks  define  the  coastline  configuration  present 
within  the  model  grid  system.  Further  output  processing 
involves  the  incorporation  of  the  initial  ocean  currents 
and  initial  surface  wind  components.  Other  input 
variables  listed  above  are  obtained  during  the  processing 
of  the  thermodynamic  module  which  accesses  the  input  mo¬ 
dule  to  obtain  the  needed  atmospheric  thermodynamic 
variables. 

The  overall  structure  of  the  processing 
model  can  also  be  considered  to  exist  in  the  three  main 
portions  defined  below.  The  first  is  a  momentum  balance 
which  includes  air  and  water  stresses,  coriolos  force, 
internal  ice  stress  and  ocean  tilt.  Non-linear  boundary 
layers  for  both  the  air-ice  and  ocean-ice  surfaces  are 
used.  A  key  component  is  the  force  due  to  internal  ice 
stress.  This  is  defined  by  a  constituitive  law  which  re¬ 
lates  the  ice  stress  to  the  strain  rate  and  ice  strength. 
For  this  model  a  viscous -plastic  constituitive  law  is 
followed. 

The  second  feature  of  the  processing 
module  consists  of  continuity  equations  describing  the 
evolution  of  the  thickness  characteristics  on  two  levels. 
Two  categories  of  ice  thickness  are  assumed;  thin  ice 
(less  than  0.5  m  in  thickness)  and  thick  ice  (greater  than 
0.5  m  in  thickness).  To  keep  account  of  these  categories 
two  variables  are  maintained;  ice  thickness  per  unit  area 
and  the  ice  concentration  which  is  defined  as  the  fraction 
of  area  of  a  grid  cell  covered  by  thick  ice.  Thermo¬ 
dynamic  terms  are  included  on  these  continuity  calcula¬ 
tions. 

The  final  component  is  an  ice  strength 
value  which  is  taken  to  depend  linearly  upon  the  ice 
thickness  and  exponentially  on  the  ice  concentration. 
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The  coupled  non-linear  equations  are 
treated  as  an  initial  value  problem  using  energy  conserv¬ 
ing  finite  differences.  The  momentum  equations  are  int¬ 
egrated  implicitly  in  order  to  avoid  a  time  limit  con¬ 
straint.  A  relaxation  technique  is  used  on  the  set  of 
simultaneous  equations  at  each  time  step. 

The  numerical  scheme  uses  a  staggered 
grid  which  allows  ice  strength  and  ice  velocities  to  vary 
in  space.  To  a  large  degree  this  staggered  grid  is 
patterned  after  those  used  in  primitive  equation  ocean 
models. 

As  mentioned  above,  initial  conditions 
at  all  points  and  ice  velocities  at  the  boundaries  are 
thereafter  required  to  initiate  the  integration  of  the 
system  of  equations.  The  most  natural  condition  is  to 
take  the  ice  velocity  to  be  zero  at  the  boundary.  This 
can  be  accomplished  at  land  boundaries  or  open  ocean 
boundaries  where  there  is  no  ice.  This  boundary  condition 
does  not  affect  the  ice  motion  in  such  circumstances 
since,  in  the  absence  of  ice  the  strength  is  zero.  It  is 
also  possible  to  set  an  "open"  boundary  condition  by  sea¬ 
ting  the  strength  equal  to  zero  near  a  boundary.  This  type 
of  open  boundary  is  used  at  the  Spitzbergen-Greenland 
passage  to  form  a  natural  inflow/outflow  region  (Figure  1) . 

In  the  computer  code,  three  boundary 
masks  are  used  to  define  the  "closed"  and  "open"  boun¬ 
daries.  Consequently  by  altering  these  masks,  highly  ir¬ 
regular  boundaries  may  be  taken  into  account. 

Because  of  the  strong  ice  interaction, 
the  momentum  equations  are  parabolic  in  form  and  hence 
have  few  numerical  instability  problems  over  longterm  in¬ 
tegrations.  To  avoid  non-linear  instabilities  in  longterm 
simulations  which  can  arise  from  the  non  linear  advection 
terms  in  the  continuity  equations,  small 
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biharmonic  and  harmonic  terms  are  added  to  the  continuity 
equation. 

A  thermodynamic  system  is  incorporated 
which  specifies  the  ice  growth  rate  as  a  function  of 
thickness  and  time  of  year. 

Input  to  ICEMDL 

Input  to  the  ICEMDL  program  exists  on  a 
number  of  files  which  must  be  set  up  before  the  actual 
model  is  run.  Only  a  portion  of  these  files  are  actually 
accessed  by  PROGRAM  ICEMDL.  The  remainder  of  the  input 
files  is  accessed  by  input  module  subroutines. 

The  input  files  are  labeled  as  TAPE7=IN 
and  TAPE8=DATE  in  the  CY203  convention.  TAPE8  contains 
the  date  time  groups  of  the  period  for  the  model  inte¬ 
gration.  The  date  time  group  ( DTG )  is  the  standard  format 
defined  at  FNOC  consisting  of: 

YYMMDDHH 

where 

YY  =  year; 

MM  =  month; 

DD  =  day; 

HH  *  hour. 


Each  DTG  is  kept  in  one  word  and  is 
valid  for  one  day  (time  step).  The  DTG  is  used  to  access 
the  proper  fields  from  the  CRANDIC  data  base  consisting  of 
the  atmospheric  data  maintained  by  FNOC.  Three  DTG  values 
are  input  to  ICEMDL  at  each  time  step  requiring  new  at¬ 
mospheric  data.  The  initial  DTG  values  input  for  the  very 
first  time  step  are  set  up  slightly  differently  than  the 
other  DTG  values.  For  all  accesses  of  the  DTG  file  the 


first  DTG  contains  the  DTG  of  the  current  day  valid  for 
that  time  step.  The  second  DTG  contains  the  DTG  for  12 
hours  before  the  current  day.  For  the  initial  group  of 
DTG 's  the  third  position  contains  the  DTG  of  the  last  time 
step  of  the  previous  run.  This  DTG  is  used  to  read  the 
final  fields  produced  by  the  previous  run.  These  fields 
are  used  to  restart  the  model  simulation.  Therefore,  if  a 
run  was  made  to  simulate  the  month  of  January  (ending 

January  31,  1981)  and  a  new  run  was  to  begin  for  February 
the  initial  line  of  the  TAPE8  file  would  be; 

81020100  81013112  81013100. 

The  DTG,  81020100,  is  the  current  DTG 

for  the  start  of  the  simulation.  The  DTG  81013100  is  the 
last  day  of  the  previous  simulation  and  will  be  used  to 
access  the  data  for  restarting  the  model.  In  reality  in¬ 
put  atmospheric  data  is  entered  into  the  model  every  4 
days  (time  steps)  .  Therefore  the  DTG  values,  held  in  the 
model,  change  in  increments  of  4  days. 

The  file  TAPE7  conatins  various  types 
of  input  data  used  by  the  model.  TAPE 7  is  accessed  by 

ICEMDL  to  set  the  following  variables; 

NX,  NY  -  dimension  of  the  grid  which  holds  the 

momentum  variables; 

NX1 , NY1  -  dimension  of  the  grid  which  holds  the 

thermodynamic  variables; 

N3  -  number  of  points  in  the  momentum 

variable  grid; 

N4  -  number  of  points  in  the  thermodynamic 

variable  grid. 


ICEMDL  Processing 


The  processing  within  ICEMDL  begins  by 
accessing  the  input  files  described  in  2. 4. 1.1.  The  next 
processing  involves  the  accessing  of  input  module  sub¬ 
routines.  Subroutine  MESH  is  the  first  input  module 
routine  which  is  accessed.  This  routine  defines  the  model 
grid  in  terms  of  the  FNOC  hemispheric  grid.  MESH  returns 
the  i,j  values  of  the  model  grid  defined  in  terms  of 
hemispheric  grid. 

Subroutine  BNDRY  is  accessed  after 
MESH.  BNDRY  returns  the  boundary  masks  to  ICEMDL. 

Subroutine  OCEAN  is  called  to  formulate 
the  u,  v  components  of  the  surface  ocean  currents  and  re¬ 
turns  them  to  ICEMDL. 

At  this  point,  some  basic  house  keeping 
tasks  are  performed  such  as  presetting  several  arrays  to 
zero  and  accessing  COPEN  which  is  an  FNOC  routine  which 
"opens"  the  CRANDIO  file  to  be  accessed  by  the  model. 

An  error  condition  develops  if  the 
routine  COPEN  finds  that  it  is  unable  to  open  the  CRANDIC 
files.  The  program  will  terminate,  stating  that  there  is 
a  COPEN  ERROR  on  the  STOP  line  of  the  dayfile. 

Subroutine  INITIAL  is  the  final  input 
module  routine  accessed  by  ICEMDL.  Two  calls  are  made  to 
obtain  the  u  and  v  atmospheric  wind  components  respective¬ 
ly. 

The  next  section  of  ICEMDL  accesses  the 
CRANDIC  data  base  written  by  previous  runs  of  the  model. 
The  variables;  ice  thickness,  ice  concentration  and  ice 
drift  are  read  to  restart  the  model  simulation. 

At  this  point  all  necessary  information 
is  available  for  the  start  of  processing.  The  i  bove  des¬ 
cribed  code  is  never  executed  again  during  the  current  run. 
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The  actual  equation  integration  starts 
with  a  series  of  calls  to  routines  FORM  and  RELAX.  FORM 
computes  all  parameters  necessary  for  use  in  the  relaxa¬ 
tion  scheme.  These  include  the  air  and  water  drag  co¬ 
efficients,  forces  terms,  ice  strength  terms,  and  the  vis¬ 
cosity  terms.  Subroutine  RELAX  performs  the  relaxation 
technique. 

There  is  a  sequence  of  two  calls  to 
these  routines.  The  first  performs  the  prediction  portion 
of  the  momentum  time  step.  Before  the  first  RELAX 
call, the  third  level  of  ice  velocities  and  the  centered 
ice  velocities,  held  in  UICE1,  VICEC.  are  set  equal  to  the 
level  of  ice  velocities  by  Hibler  (1980) .  During  this 
procedure,  the  time  step  is  halved.  FORM  is  called  in 
this  first  step  to  use  the  present  ice  velocity  values  to 
linearize  the  momentum  equations. 

The  second  calls  of  FCRM  and  RELAX 
amount  to  the  main  forward  time  step  of  the  "corrector" 
section  of  this  "predictor-corrector"  method.  In  this 
case  the  "predicted"  values  of  the  ice  velocities  are  used 
in  the  second  FORM  call  to  estimate  the  viscosity  para¬ 
meters. 

After  the  momentum  equations  have  been 
implicitly  stepped  foreward  using  the  relaxation  tech¬ 
nique,  several  diagnostic  calculations  are  carried  out. 
The  squared  ice  velocity  and  the  squared  ice  velocity 
difference  between  the  times,  t  and  t+1  are  computed. 
This  is  a  simple  measure  of  the  change  taking  place  in  the 
ice  velocity  field  during  the  time  step  advance. 

The  predicted  ice  velocity  values  are 
contained  in  the  first  level  of  the  UICE  and  VICE  arrays. 

The  ice  velocity  and  divergence  values 
are  passed  to  the  output  module  routines,  UVPLOT  and 
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DIVERG  respectively. 

Following  this  momentum  time  step,  the 
thermodynamic  equations  are  explicitly  stepped  forward  in 
time.  Subroutine  ADVECT  is  called  to  handle  the  dynamical 
portions  of  the  continuity  equations  for  ice  thickness  and 
concentration.  Subroutine  HEAT  is  called  to  obtain  the 
thick  and  thin  ice  growth  rates  through  the  use  of  a  heat 
budget.  Subroutine  GROWTH  is  called  to  heat  the  thermo¬ 
dynamic  portion  of  the  continuity  equations. 

This  concludes  processing  for  the  time 
integration.  the  remainder  of  the  code  is  used  to  monitor 
diagnostic  variables  and  the  output  modules. 

Subroutine  XSUM  is  used  to  obtain  the 
total  ice  held  within  the  grid,  excluding  outflow  cells. 
The  following  diagnostic  values  are  then  computed; 

i)  Total  open  water  growth  for  each 
grid  cell;  HDIFF; 

ii)  Net  open  water  growth  for  the 
basin;  GRSUM1; 

iii)  Net  ice  growth  for  the  basin; 

FHS UM ; 

iv)  Total  ice  in  the  outflow  cells; 
TOUT. 

The  ice  held  in  the  outflow  cells  is 
explicitly  determined  through  variables,  THEFF  and 
THEFF1 .  The  variable  THEFF1  contains  the  amount  of  ice  in 
the  open  cells  and  is  computed  at  the  beginning  of  the 
time  step. 

The  output  module  begins  with  a  list  of 
PRINT  statements  which  form  the  hard  copy  output.  Hard 
copy  output  is  printed  every  four  time  steps.  The 
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variable,  LSTEP,  is  used  to  count  time  steps  and  branch  to 
the  output  section  on  the  fourth  time  step. 

Subroutine  PRNT  is  used  to  print  the 

model  arrays. 


Subroutine  ADJUST  is  called  after  the 
printed  output  is  formulated.  ADJUST  is  used  to  define  the 
amount  of  ice  held  in  the  outflow  cells  for  use  at  the 
beginning  of  the  next  time  step. 

Subroutine  GROWDEC  and  HAPLOT  are  used 
in  the  output  module  to  output  the  ice  thickness,  ice  con¬ 
centrations  and  growth/decay  rates  to  the  CRANDIO  file 
maintained  by  the  model. 

The  final  function  performed  in  ICEMDL 
is  to  decide  if  processing  is  complete.  The  variable 
ITSTEP  which  is  input  from  TAPE8  defines  the  total  number 
of  time  steps  minus  one  to  be  used  in  the  run.  The 
variable  ICOUNT  is  used  to  keep  track  of  how  many  steps 
have  currently  been  executed. 

If  more  time  steps  are  required,  a 
check  is  made  to  determine  whether  new  atmospheric  wind 
data  is  needed.  New  wind  data  is  used  every  four  days, 
the  variable,  LSTEP,  defines  the  fourth  time  step  as  des¬ 
cribed  above.  Subroutine  INITIAL  is  called  to  access  the 
CRANDIO  file  containing  the  wind  data  if  needed.  The  pro¬ 
gram  continues  processing  until  time  steps  are  no  longer 
desired . 


If  the  number  of  desired  time  steps  has 
been  completed,  the  diagnostic  variables  of; 


i)  outflow; 

ii)  net  ice  growth; 

iii)  net  open  water  growth; 


are  written  to  the  file  TAPE3  for  use  in  a  restart  run  if 
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desired.  Subroutine  STATPRT  is  called  to  print  time  use 
statistics  on  various  routines  used. 


ICEMDL  Output 

The  output  of  ICEMDL  consists  of  fore¬ 
casts  of  the  following  variables; 

i)  ice  thickness; 

ii)  ice  concentration; 

iii)  ice  drift; 

iv)  ice  divergence/convergence; 

v)  ice  strength; 

vi)  ice  growth. 

These  variables  are  output  in  printed 

form  and  also  in  CRANDIC  files  on  the  Cy203.  CRANDIG  is 

the  type  of  file  used  operationally  at  FNOC,  on  the  Cy203, 
for  maintenance  of  the  environmental  data  base.  CRANDIC 
is  analoguous  to  ZRANDIG  on  the  Cyl70's  and  6600's  at  FNOC. 

The  output  CRANDIO  files  contain  one 
record,  produced  every  four  time  steps,  for  each  of  the 
above  6  variables.  Specific  information  as  to  the  format 
and  structure  of  a  CRANDIG  file  can  be  found  in  the 

appropriate  FNOC  technical  write-up. 

The  records  are  labeled  with  the  cata¬ 
log  name  of  each  variable.  The  date  of  the  record,  and  a 
tau  value.  The  catalog  names  are  defined  as; 


i) 

FFF  - 

ice 

ii) 

DDD  - 

ice 

iii) 

THK  - 

ice 

iv) 

CON  - 

ice 
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speed ; 
direction; 
thickness ; 
concentration; 


i 


1 


v)  PRS  -  ice  strength; 

vi)  HDF  -  ice  growth; 

vii)  GAR  -  open  water  growth; 

viii)  DIV  -  convergence/divergence. 

Interfaces 

Program  ICEMDL  interfaces  with  sub¬ 
routines  which  comprise  the  input,  processing  and  output 
modules.  The  following  table  defines  all  interfaces  bet¬ 
ween  these  routines  and  ICEMDL.  Specific  arrays  are  de¬ 
fined  under  Tables  and  Items. 

Subroutine  Name  Interfaces 

MESH  receives  -  NX,  NY,  NXl,  NY1 

NUMBER  -  number  of  points  in 
grid; 

returns  -  GRDI  -  i  grid  points; 

GRDJ  -  j  grid  points; 

BNDRY  receives  -  NX,  NY,  NXl,  NY1 

returns  HEFFM  -  thermodynamic 

variables  boundary  mask; 

UVM  -  momentum  variables 
boundary  mask; 

OUT  -  outflow  boundary  mask; 
OCEAN  receives  -  NX,  NY,  GRDI,  GRDJ,  UVM; 

returns  -  GWATX  -  ocean  current  u 
components 

GWATY  -  ocean  current  v 
components ; 

COPEN  receives  -  IFILE  -  CRANDIO  file  name 

returns  -  ISTAT  -  status  of  file  open 
attempt ; 
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INITIAL 

receives 

file  unit  number,  NUMBER, 

GRDI,  GRDJ ,  IDTG ,  DTG  array. 

ITAU  -  tau  value; 

returns 

variables  read  from  FNOC  dat< 

base ; 

C READER 

receives 

IFILE  -  array  for  data; 

LABEL  -  record  name; 

catalog  name; 

record  length; 

returns 

read  status,  IS; 

DDFFUN 

receives 

DD,  FF  -  direction  and  force 

of  current  or  drift; 

returns 

U,  V  components; 

XSUM 

receives 

HEFF  -  ice  thickness  array 

NX1 ,  NY 1 ; 

returns 

-  total  of  array  -  THEFF  or 

THEFF1 ; 

PRNT 

receives 

array  name; 

dimensions  of  array; 
positions  to  be  printed; 

FORM 

receives 

-  UICE ,  VICE,  ETA,  ZETA,  AMASS, 

GAIRX ,  GAIRY ,  GWATX,  GWATY, 

OUT,  HEFFM ,  NX,  NY,  NXl,  NY1 , 

HEFF,  AREA; 

returns 

-  DRAGS,  DRAGA,  DIV; 

RELAX 

receives 

-  UICE,  VICE,  ETA,  ZETA,  AMASS, 

GAIRX,  GAIRY,  GWATX,  GWATY, 

DRAGS,  DRAGA,  OUT,  HEFFM,  NX, 

NY,  NXl,  NY1,  HEFF,  AREA; 

returns 

-  UICE,  VICE; 

DIVERG 

receives 

-  DIV,  NX,  NY,  GRDI,  GRDJ,  IDTG, 

ITAU; 

UVPLOT 

rece ives 

-  UICEC ,  VICEC,  GRDI,  GRDJ,  NX, 

NY; 
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ADVECT 

receives 

-  NICEC,  VICEC ,  HEFF,  DIFFI, 
LAD,  HEFFM ,  NX,  NY,  NX1,  NYl ; 

returns 

-  HEFF  or  AREA,  DIFFI; 

HEAT 

receives 

-  GRDI,  GRDJ ,  HEFF,  AREA,  GAIRX 

GAIRY,  ITAU,  IDTG ,  NXl ,  NYl, 

NUMBER; 

returns 

-  FC,  FHEFF ; 

GROWTH 

receives 

-  HEFF,  AREA,  HC,  A22,  FHEFF, 

FO,  HCORR ,  HEFFM,  OUT,  NXl, 
NYl; 

returns 

-  HEFF,  AREA,  GAREA; 

HAPLOT 

receives 

-  HEFF,  AREA,  IDTG,  ITAU,  GRDI, 

GRDJ,  NX,  NY; 

GRCKDEC 

receives 

-  HDIFF ,  FHEFF,  GAREA,  IDTG, 

ITAU,  NXl,  NYl; 

ADJUST 

receives 

-  HEFF,  AREA,  OUT,  HEFFM,  NX, 

NY,  NXl,  NYl; 

returns 

-  HEFF,  AREA. 

ICEMDL  Tables  and  Items 

The  following  lists  define  all  common 
blocks  and  major  variable  items  used  in  ICEMDL. 

I.  Common  Blocks 

/BUOY/  BUOYI,  BUCYJ  -  Buoy  position  in  terms  of  the  FNGC 

I,  J  grid; 

BX,  BY  -  Buoy  positions  in  terms  of  the 

model  x,  y  grid; 

/FORCE/  FORCEX  -  x  component  of  external  force  plus 

ice  pressure  gradient; 
y  component  of  external  force  plus 
ice  pressure  gradient; 
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/STEP/ 

DELTAT 

Time  step  in  seconds; 

DELTAX, 

mesh  size  in  meters,  x,  y 

DELTAY 

directions  respectively; 

/PRESS/ 

IDENT 

CRANDIC  record  indentification 

block ; 

DATA 

hemispheric  atmospheric  data; 

FILL 

filler  to  put  the  block  on  small 
page  boundary  (required  by  CRANDIC 

software) ; 

/IJ/ 

IJ63 

I,J  grid  points  of  SKILES  current 

data ; 

/DD/ 

ID 

CRANDIC  identification  block  for 

ice  drift  direction; 

DD 

ice  drift  directions; 

FILLD 

filler  for  small  page  boundary; 

/FF/ 

w 

1 

Same  as  /DD/  except  for  ice  drift 

FILLF 

speeds; 

/AR/ 

IDA,  ARRAY,  - 

Same  as  /DD/  except  for  ice 

FL 

concentration. 

Variable  Items  : 


UICE 

- 

u  component  of  ice  velocity 

VICE 

- 

v  component  of  ice  velocity 

ETA 

- 

non  linear  shear  viscosity 

ZETA 

- 

non  linear  bulk  viscosity 

GAIRX 

- 

u  component  of  wind 

GAIRY 

- 

v  componet  of  wind 

AMASS 

- 

ice  mass  per  unit  area 

HEFF 

- 

mean  ice  thickness  per  unit  area 

AREA 

- 

fraction  of  area  covered  by  thick 

UICEC 

- 

Intermediate  ice  velocities  for 

VICEC 

- 

use  in  semi-implicit  time  step 

GWATX 

- 

u  component  of  ocean  currents 

GWATY 

- 

v  component  of  ocean  currents 
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STRESS 

FHEFF 

FO 

DRAG  A 
DRAGS 
DIV 


XX,  YY  and  XY  components  of  ice  stress 

growth  rate  of  thick  ice 

growth  rate  of  thin  ice 

water  drag  plus  Coriolos  parameter 

water  drag  plus  inertial  term 

ice  convergence/divergence 


2.4.2 


Input  Module  Subroutines 


2.4. 2.1  Subroutine  MESH 

Subroutine  MESH  is  used  to  create  the 
model  grid.  The  grid  is  constructed  as  a  subset  of  the 
FNCC  hemispheric  grid. 


Input  to  MESH 

Subroutine  MESH  receives  input  from  2 
sources.  The  first  source  is  the  formal  parameters  from 
ICEMDL.  MESH  also  reads  the  input  file  labeled  TAPE7 . 
The  following  variables  are  contained  in  the  formal  para¬ 
meter  list; 


i)  GRDI,  GRDJ  -  computed  in  MESH, 

and  contain  the  i, 
j  coordinates  of  the 
model  grid 

ii)  NX,  NY,  NX1,  NYl ,  N  -  grid  sizes, 

defined  in  TAPE  7. 

The  following  variables  are  obtained 
from  the  input  file  TAPE7; 
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i) 

10, 

11  - 

Defining  i  grid  points  on 
the  FNOC  hemispheric  grid; 

ii) 

JO, 

J1  - 

Defining  j  grid  points  on 
the  FNOC  hemispheric  grid; 

iii) 

N 

The  number  of  points  in 

each  row  and  column  in  the 

model  grid. 

Subroutine  MESH  creates  a  square  grid. 
The  variables  10,  II,  JO,  J1  are  defined  as  follows; 

x  x 

i,  j  II#  J1 

x  x 

10, JO  i ,  j 

Therefore  10,  JO  define  the  bottom  left 
corner  of  the  desired  model  grid  and  II,  J1  define  the 
upper  right  corner  of  the  model  grid.  MESH  fills  in  the 
remainder  of  the  grid,  depending  upon  N.  All  i,  j  coor¬ 
dinates  are  in  reference  to  the  hemispheric,  63  x  63  grid 
of  FNOC. 


Processing 

The  processing  of  subroutine  MESH  be¬ 
gins  by  reading  all  necessary  input  data  from  TAPE7 .  From 
the  variable,  N,  the  values  of  NX,  NY  and  NXl,  NY1  are 
computed.  The  actual  grid  is  then  formulated  using  the 
specified  corner  points  and  the  number  of  points  desired 
in  each  row  (column).  The  column  points  are  defined, 
stored  in  GRDJ,  followed  by  the  row  points  stored  in  GRDI. 
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After  the  entire  grid  is  formed,  the 
MESH  size  is  computed  utilizing  the  characteristics  and 
map  factor  of  the  polar  stereographic  grid.  The  map  fac¬ 
tor  of  the  model  grid  is  also  computed. 

The  final  function  of  the  routine  is  to 
compute  the  initial  buoy  positions  in  terms  of  the  newly 
constructed  model  grid. 

Output 

Subroutine  MESH  produces  printed  output 
specifying  the  following; 

i)  Corner  grid  points; 

ii)  mesh  size  of  model  grid; 

iii)  map  factor; 

iv)  buoy  grid  locations. 

Interfaces 

Subroutine  MESH  does  not  call  any  other 

subroutines. 

Tables  and  Terms 

All  common  blocks  and  major  variables, 
used  in  MESH  are  defined  under  the  ICEMDL  section. 

2.4. 2. 2  Subroutine  BNDRY 


Subroutine  BNDRY  is  called  to  set  up 
the  boundary  masks  for  the  thermodynamic,  momentum  and 
outflow  grids.  A  boundary  mask  consists  of  an  array  which 
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contains  a  1  in  grid  locations  where  computations  are  per¬ 
formed  and  a  0  on  all  boundary  points.  By  altering  these 
masks  one  can  obtain  any  desired  boundary  or  coastline 
configuration. 


Input 

The  respective  boundary  masks  are  read 
from  the  input  file  labeled  TAPE7.  The  grid  sizes  are 
input  through  the  formal  parameter  list  interfaced  with 
ICEMDL. 


Processing 

The  standard  FORTRAN  READ  function  is 
used  to  access  TAPE7  and  move  the  boundary  masks  to  the 
respective  locations. 

Output 

The  arrays,  UVM,  HEFFM  and  OUT  are  cre¬ 
ated  by  BNDRY . 


Interfaces 


Subroutine  BNDRY  does  not  call  any 

other  subroutines. 


Tables  and  Items 

The  following  major  variables  are  used; 

UVM  -  boundary  mask  for  momentum 
variables ; 
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HEFFM 


boundary  mask  for 
thermodynamic  variables 
OUT  -  boundary  mask  containing 
outflow  points. 


2. 4. 2. 3  Subroutine  OCEAN 

Subroutine  OCEAN  is  used  to  define  the 
surface  ocean  currents.  These  currents  are  defined  by  the 
SKILES  sea  ice  drift  model  used  at  FNOC. 

Input 

Ocean  current  directions  are  read  from 
the  input  file  TAPE7 .  These  direction  values  are  on  the 
model  grid.  The  values  were  interpolated  from  the  SKILES 
grid  to  the  current  grid  and  the  results  placed  on  TAPE7. 

The  model  grid  size  and  i,  j  coor¬ 

dinates  of  the  model  grid  are  passed  to  OCEAN  from  ICEMDL 
through  the  formal  parameter  list. 

Processing 

The  ocean  current  magnitudes  are  lo¬ 
cated  in  the  DATA  statement  defining  the  array  WF.  The 

ocean  current  magnitudes  are  input  from  TAPE7  through  the 

use  of  a  standard  READ.  The  current  magnitudes  are  con¬ 
verted  from  cm/sec  to  m/sec  and  placed  into  a  model  grid. 
The  current  directions,  on  TAPE7 ,  also  needed  to  be  multi¬ 
plied  by  10.  For  example  a  direction  of  270  is  read  in  as 
27. 

The  current  directions  and  magnitudes 
are  converted  to  u,  v  components  with  the  FNOC  routine 
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DDFFUV .  Details  pertaining  to  the  function  of  this 
routine  are  contained  in  the  standard  FNOC  subroutine 
write-ups.  The  code  is  placed  within  this  program  because 
the  library  is  not  available  to  the  CY203  at  the  time  of 
this  writing. 

Output 


The  arrays  GNATX,  GVvATY  are  produced  in 

OCEAN. 

Interfaces 

Subroutine  OCEAN  interfaces  with  sub¬ 
routine  DDFFlw  which  is  a  standard  FNOC  library  routine 
used  to  convert  a  direction  and  magnitude  to  u,  v  com¬ 
ponents  . 

Tables  and  items 


All  major  variables  contained  in  OCEAN 
are  defined  under  ICEMDL. 


2. 4. 2. 4 


Subroutine  INITIAL 


Subroutine  INITIAL  is  used  to  access 
FNOC  CRANDIC  data  on  the  CY203  and  place  the  specified 
atmospheric  data  into  the  model  grid. 


Input 


The  main  input  to  INITIAL  is  obtained 
from  ICEMDL  through  the  formal  parameter  list.  The 


variable,  NO,  specifies  which  field  is  to  be  accessed  from 
the  CRANDIO  data  base.  Value  of  NO  specifies  which  posi¬ 
tion  in  array  IRCD  is  to  be  used.  IRCD  holds  the  catalog 
names  of  the  various  atmospheric  data  input  to  the  model. 
The  array,  IDTG,  contains  the  data  time  group  to  be  used 
in  reading  the  data.  ITAU  is  the  tau  value  and  N  is  the 
number  of  rows  (columns)  in  the  model  grid.  The  arrays 
GRDI,  GRDJ  are  defined  as  in  MESH. 

Processing 

The  first  section  of  code  in  INITIAL  is 
setting  up  the  CRANDIG  record  name.  This  is  a  2  word 
ASCII  label.  The  first  word  contains  the  following; 

i)  Catalog  name  of  data; 

ii)  flaps  character; 

iii)  length  of  record. 

The  second  word  contains  the  DTG.  All 
of  this  information  is  masked  together  in  array  LABEL. 
The  FNCC  routine,  CHECKNC,  is  used  to  determine  whether 
the  specified  data  is  present  on  the  CRANDIG  file.  If  not 
present  the  program  terminates  with  the  following  message; 

STOP  CHECKNC  NO  DATA  INITIAL. 

Providing  the  data  is  present,  the  FNOC 
routine  CREADER  is  used  to  read  the  hemispheric  data  into 
array,  DATA. 

Subroutine  INTRP  is  used  to  interpolate 
the  hemispheric  data  to  the  model  grid.  INTRP  is  an  FNOC 
subroutine  which  performs  the  interpolating  function. 
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A  detailed  description  of  the  operations  of  INTRP  is 
available  in  the  FNOC  utility  subroutine  documentation. 

Outputs 

Subroutine  INITIAL  provides  for  a  hard 
copy  print  of  the  requested  record  label.  The  array, 
GAIR,  is  used  to  store  the  resultant  data  placed  in  the 
model  grid. 


Interfaces 


Subroutine  INITIAL  interfaces  with  the 
CRANDIO  data  base  software  at  FNOC  plus  a  utility  library 
routine,  INTRP.  Detailed  descriptions  on  all  these 
routines  are  available  as  standard  products  of  FNOC. 

Tables  and  Items 


The  major  variables,  used  in  INITIAL, 
are  defined  in  ICEMDL.  The  array  IRCD  is  used  to  maintain 
the  respective  catalog  names  of  the  CRANDIO  data.  These 
catalog  names  are  defined  as  follows; 


i) 

A20 

u  wind  components; 

ii) 

A21 

v  wind  componments; 

iii) 

A10 

atmospheric  temperature  at 
the  surface; 

iv) 

A01 

surface  pressure; 

v) 

A12 

surface  vapor  pressure; 

vi) 

All 

shortwave  radiation; 

vii) 

A18 

-  total  heat  flux; 

viii) 

A16 

sensible  plus  evaporative 
heat  flux. 
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2.4.3 


2. 4. 3. 2  Subroutine  FORM 

Subroutine  FORM  is  used  to  calculate 
the  drag  coefficients,  external  forces  and  ice  strength 
parameters  for  use  in  the  time  integration  of  the  momentum 
equations . 
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Input 


Subroutine  FORM  receives  all  required 
input  from  ICEMDL  through  the  formal  parameters.  All  of 
the  variable  definitions  passed  to  FORM  are  presented  un¬ 
der  ICEMDL.  Certain  constants  are  defined  as  follows; 


i) 

FCOR 

-  average  coriolis 

parameter ; 

ii) 

RHGAIR 

-  air  density; 

iii) 

SINWIN 

-  sine  of  air  turning 
angle ; 

iv) 

COSWIN 

-  cosine  of  air  turning 
angle ; 

V) 

SINWAT 

-  sine  of  ocean  turning 
angle  ? 

vi) 

COSWAT 

-  cosine  of  ocean  turning 
angle; 

vii) 

ECCEN 

-  ratio  of  the  axes  of 

the  plastic  yield 
curve. 


Processing 

Subroutine  FORM  initially  calculates 
the  variables  required  to  obtain  the  external  forces  oper¬ 
ating  upon  the  ice.  Within  the  first  major  loop  the  ice 
mass  per  unit  arrea,  coriolis,  non-linear  water  stress 
coefficient  and  antisymmetric  water  drag,  DRAGA,  are  com¬ 
puted.  The  next  major  loop  calculates  the  non-linear  air 
stress  and  the  symmetric  water  drag  term,  DRAGS.  These 
terms  are  calculated  in  a  separate  loop  due  to  a  different 
grid  requirements  of  the  GAIRX,  GAIRY  (wind  component) 
term  (See  Appendix  A) . 
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The  remaining  portion  of  the  subroutine 
deals  with  finalyzing  the  external  force  terms  and  access¬ 
ing  Subroutine  PLAST  to  compute  non-linear  shear  and  bulk 
viscosities  respectively.  Before  PLAST  is  called,  the  ice 
strength,  PRESS,  is  computed.  After  PLAST  is  called  the 
computed  viscosities  are  set  to  zero  at  the  outflow  points 
by  multiplying  by  the  boundary  mask  OUT. 

Finally,  the  external  force  components 
term  is  combined  with  the  ice  pressure  gradient. 


The  results  of  processing  in  routine 
FORM  are  stored  in  DRAGS,  DRAGA,  FGRCEX ,  FORCEY,  ETA, 
ZETA,  and  PRESS.  These  variables  are  defined  under  ICEMDL. 

Interfaces 

Subroutine  FORM  interfaces  with  routine 
PLAST.  Subroutine  PLAST  calculates  the  non  linear  vis¬ 
cosities  based  on  a  Plastic  yield  curve. 

Tables  and  Items 

Subroutine  FORM  operates  upon  a  number 
of  major  variables  which  are  defined  under  ICEMDL. 

2. 4. 3. 3  Subroutine  PLAST 

Subroutine  PLAST  is  used  by  subroutine 
FORM  to  calculate  the  non  linear  viscosity  terms  based  on 
a  plastic  yield  curve. 


1 


I  nput 


Subroutine 
parameters  from  FORM,  through 


PLAST  receives  all  input 
the  formal  parameter  list. 


Processing 

The  first  main  loop  of  PLAST  uses  the 
ice  drift  components  to  calculate  the  XX,  XY,  YY  strain 
rates  of  the  ice  (Ell,  E12,  E22  respectively).  These  com¬ 
ponents  are  used,  with  the  constituitive  law  to  calculate 
the  non-linear  bulk  viscosity.  The  bulk  viscosity  is  in 
turn  used  to  compute  the  non-linear  shear  viscosity. 

The  final  computation  within  PLAST  is 
the  calculation  of  the  XX,  XY,  YY  components  of  ice  stress 
plus  the  divergence  is  calculated  from  the  XX,  YY  strain 
rates . 

Output 


The  main  output 
in  the  arrays  ETA,  ZETA  and  sent  to 
meter  list. 


of  PLAST  is 
FORM  through 


contained 
the  para- 


Inter faces 


routines. 


Subroutine  PLAST  calls  no  other  sub¬ 


Tables  and  Items 


The  main  variables  of  PLAST  which  have 
not  been  identified  previously  are; 

i)  Ell,  E12,  E22  -  XX,  XY,  YY  strain 

rates  respectively. 
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2. 4. 3. 4 


Subroutine  RELAX 


Subroutine  RELAX  is  the  main  routine  of 
the  processing  module.  This  routine  applies  a  relaxation 
technique  to  the  dynamical  equations  for  their  numerical 
integration  in  time.  Much  numerical  detail  is  contained 
in  the  routine  and  described  by  Hibler  (1979)  . 

Input 

Subroutine  RELAX  receives  all  input 
from  ICEMDL  through  the  parameter  list. 

Processing 

The  processing  of  RELAX  is  broken  into 

a  number  of  separate  modules.  The  first  3  major  loops 

perform  operations  involving  the  previous  time  values  of 
u,  v  plus  the  evaluation  of  the  diagonal  components  of  the 
computation  matrix. 

The  main  loop  (103)  performs  the  iter¬ 
ative  relaxation  scheme.  This  loop  performs  all  necessary 
calculations  to  obtain  a  new  estimate  of  the  u,  v  ice 
drift  components  (UICE,  VICE) .  After  the  new  components 
are  calculated,  2  checks  are  made.  The  first  check  ex¬ 
amines  the  number  of  iterations  completed  to  determine  if 
more  than  1300  have  taken  place.  If  so,  the  routine  shall 
end  printing  a  message  stating  more  than  1300  iterations 

have  occurred  with  no  convergence.  The  second  check 
searches  for  the  100th  iteration.  At  this  point  the 
routine  switches  from  a  over  relaxation  scheme  to  a 
straight  relaxation  scheme. 

After  the  checks  have  been  executed, 

the  difference  between  the  new  solution  and  previous 
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iterative  solution  is  computed.  if  the  difference  lies 
within  an  accepted  tolerance  the  routine  ends.  If  the 
difference  does  not  meet  the  tolerance  specification 
another  iteration  is  performed.  The  old  iteration  value 
is  stored  in  the  third  level  of  arrays  UICE,  VICE  while 
the  new  iterative  solution  is  placed  in  the  first  level  of 
arrays. 

The  relaxation  code  is  made  more  com¬ 
plex  by  the  separation  of  the  finite  difference  computa¬ 
tions  into  subroutines,  FELLIP  and  FELLD1 .  The  code  is 
also  generalized  to  fit  into  the  predictor-corrector 
scheme  previously  defined.  The  parameter,  THETA,  defines 
whether  backwards,  centered,  or  forward  time  steps  are 
used.  A  value  of  0.5  initiates  a  centered  time  and  a 
value  of  0  dictates  a  forward  step. 

Outputs 

A  printed  output  message  is  made  at  the 
completion  of  the  relaxation  procedure.  The  message 
states  the  number  of  iterations  used  and  the  value  of  the 
difference  between  iterations  at  the  end.  The  results  of 
the  processing  are  placed  in  UICE,  VICE,  level  1. 

Interfaces 

Subroutine  RELAX  interfaces  with 
routines  FELLIP  and  FELLD1.  These  routines  calculate 
finite  differences  used  in  the  relaxation  routine. 


Tables  and  Items 

A  large  number  of  variables  are  used 
internally  by  RELAX.  All  major  variables,  however,  are 
defined  under  the  description  of  ICEMDL. 
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Subroutine  FELLIP 


2. 4. 3. 5 


Subroutine  FELLIP  is  used  to  calculate 
finite  differences  used  in  the  ralaxation  technique. 

Input 

All  required  input  to  FELLIP  is  pro¬ 
vided  in  the  parameter  list  and  passed  to  FELLIP  from 
RELAX. 


Processing 

Subroutine  FELLIP  operates  on  one  grid 
position  at  a  time.  The  position  is  defined  by  the  input 
variable,  i,  j,  k.  FELLIP  is  called  a  number  of  times 
producing  various  terms  needed  by  the  relaxation  routine 
at  each  call. 

The  code  remains  a  very  straight  for¬ 
ward  calculation  of  the  finite  difference  approximations 
of  the  specified  terms. 


Output 

The  array,  F,  holds  all  resultant  cal¬ 
culations. 


Interfaces 

No  further  subroutines  are  called  by 

FELLIP. 

Tables  and  Items 

No  new  major  variables  are  used  in 


FELLIP. 
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2. 4. 3. 6 


Subroutine  FELLD1 


Subroutine  FELLD1  is  also  used  to  cal¬ 
culate  finite  differences  for  the  relaxation  code. 

Inputs 

All  required  input  to  FELLD1  is 
supplied  in  the  parameter  list  by  RELAX. 


Processing 

.  Unlike  FELLIP,  FELLD1  operates  upon  the 
entire  grid  during  one  call.  The  level  of  the  array  to  be 
altered  is  specified  by  input  variable,  K  (e.g.,  1,  2,  or 
3).  The  input  variable,  SI,  defines  the  differencing  con¬ 
stants. 


Outputs 

The  output  of  FELLD1  is  held  within  the 

array,  F. 

Interfaces 

No  subroutines  are  called  by  FELLD1. 
Tables  and  Items 

No  new  major  variables  are  used  in 

FELLD1. 
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2.4. 3. 7 


Subroutine  ADVECT 


Subroutine  ADVECT  handles  the  dynamical 
portions  of  the  thermodynamic  continuity  equations. 

Input 

Subroutine  ADVECT  receives  all  ne¬ 
cessary  input  from  ICEMDL  through  the  formal  parameter 
list.  These  variables  are  defined  under  ICEMDL  and  any 
functions  performed  by  these  variables  are  defined  below. 

Processing 

Subroutine  ADVECT  performs  the  explicit 
time  stepping  procedure  of  the  dynamical  portions  of  the 
thermodynamic  continuity  equations.  The  routine  is  de¬ 
signed  to  operate  in  two  separate  finite  difference 
schemes.  The  input  variable,  LAD,  determines  whether  a 
backward  Euler  or  Leapfrog  scheme  is  followed.  If  LAD  is 
1.0  then  the  leapfrog  scheme  is  followed,  otherwise  the 
Backward  Euler  is  used. 

Initially  the  ice  thickness  (HEFF)  and 
ice  concentrations  (AREA)  are  stepped  forward  in  time  by 
transferring  the  grid  point  values  to  the  next  lower 
level.  Therefore,  the  current  values  are  moved  to  level  2 
of  the  array  while  the  new  values  are  put  into  level  1. 
The  subroutine  is  called  twice  by  ICEMDL,  one  time  pro¬ 
cessing  HEFF  and  secondly  processing  AREA. 

The  finite  difference  approximation  to 
the  respective  variable  and  a  following  check  on  the 
finite  differencing  scheme  are  the  next  major  processing 
actions.  If  the  Backward  Euler  scheme  is  used  (LAD  is  0) 
the  scheme  is  continued  to  finish  the  necessary 
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computations  required  by  this  scheme. 

After  the  computation  is  complete  for 
both  shemes,  the  subroutine  DIFFUS  is  called  to  calculate 
a  smoothing  term  which  is  used  to  comprise  the  final  data 
value. 


The  array  HEFF,  which  is  internal  to 
ADVECT,  holds  the  final  calculations.  These  values  are 
trasferred  to  HEFF  and  AREA  in  ICEMDL. 

Interfaces 

Subroutine  ADVECT  utilizes  routine 
DIFFUS  for  calculation  of  a  smoothing  operator  used  to 
reduce  the  effort  of  non  linear  instabilities  arising  from 
non  linear  terms  in  the  continuity  equations. 

Tables  and  Items 

No  major  variables  are  used  in  ADVECT 
which  have  not  been  previously  defined. 

2. 3. 3. 8  Subroutine  DIFFUS 

Subroutine  DIFFUS  is  used  to  calculate 
small  diffusion  terms  which  are  used  to  reduce  in¬ 
stabilities  within  the  non  linear  continuity  equation. 


Inputs 

All  required  input  for  DIFFUS  is  input 
to  the  routine  by  ADVECT  through  the  formal  parameter 
list.  Major  input  variables  are  defined  under  ICEMDL. 
Any  important  functions  performed  by  these  variables  in 
DIFFUS  are  detailed  below. 

Processing 

Subroutine  DIFFUS  applies  a  simple 
smoothing  operator  to  obtain  a  small  diffusion  term  for 
the  respective  field  being  analyzed.  This  term  is  applied 
to  the  results  of  the  computations  of  the  dynamical  por¬ 
tions  of  the  continuity  equations. 

Output 

The  diffusion  terms  are  stored  in  the 
third  level  of  the  array,  HEFF. 

Interfaces 


DIFFUS . 


No  major  subroutines  are  accessed  by 


Tables  and  Items 


Subroutine  DIFFUS  introduces  no  pre¬ 
viously  defined  variables. 
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2.4. 3. 9 


Subroutine  HEAT 


Subroutine  HEAT  is  the  driving  routine 
for  the  heat  budget  code.  This  budget  solves  for  the 
thermodynamic  growth  rate  of  the  thick  and  thin  ice. 

Input 

Subroutine  HEAT  receives  all  input 
variables  from  ICEMDL,  through  the  formal  parameter  list. 

Subroutine  HEAT  also  accesses  the 
CRANDIO  data  base  which  contains  the  FNOC  environmental 
data.  The  following  data  records  are  accessed; 

i)  air  temperature; 

ii)  surface  pressure; 

iii)  surface  vapor  pressure; 

iv)  incoming  solar  radiation; 

v)  sensible  heat  flux; 

vi)  sensible  plus  evaporative  heat 
flux. 

The  data  is  read  from  the  CRANDIO  file 
through  the  use  of  subroutine  INITIAL  of  the  input  module. 

Processing 

The  primary  function  of  subroutine  HEAT 
is  to  set  up  all  necessary  variables  for  the  heat  budget 
calculations,  performed  in  subroutine  BUDGET. 

The  wind  data  is  calculated  in  the 
first  main  loop  of  HEAT.  This  variable  is  calculated  from 
GAIRX  and  GAIRY  which  contain  the  u  and  v  components 
respectively. 
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Subroutine  INITIAL  is  called  to  read 
the  air  temperature  data  at  the  surface  from  the  CRANDIO 
data  base.  Subroutine  AVG  is  called  to  compute  the  grid 
cell  average  of  the  air  temperature  data.  Subroutine  AVG 
is  used  for  all  thermodynamic  variables,  within  HEAT  to 
define  them  within  the  staggered  grid.  The  temperature 
data  is  finally  converted  from  centigrade  to  Kelvin. 

Subroutine  INITIAL  and  AVG  are  again 
used  to  define  the  atmospheric  surface  pressure  and  vapor 
pressure.  These  variables  are  used  to  calculate  the 
moisture  at  the  surface  which  is  stored  within  QA. 

Subroutine  INITIAL  is  finally  used  to 
retrieve  the  shortwave  radiation,  sensible  heat  flux  and 
sensible  plus  evaporative  heat  flux.  The  variables  are 
converted  from  the  CGS  system  to  the  MKS  system  and  used 
to  calculate  the  incoming  longwave  radiation.  At  the  end 
of  the  processing  the  following  variables  have  been  set  up 
for  use  in  BUDGET; 


i) 

TAIR 

-  air  temperature; 

ii) 

QA 

-  surface  moisture; 

iii) 

FSH 

-  incoming  shortwave 

radiation ; 

iv) 

FLC 

-  incoming  longwave 

radiation. 

The  final  preparation  before  BUDGET  is 
called  is  the  definition  of  the  mixed  layer  depth  (HMIX) 
plus  the  temperature  of  the  mixed  layer  and  temperature  of 
the  ice,  TMIX  and  TICE  respectively. 

The  variable,  KOPEN,  is  used  as  a  flag 
to  BUDGET  to  determine  whether  thin  ice  or  thick  ice 
growth  rates  are  evaluated.  Subroutine  BUDGET  is  then 
called  to  calculate  the  growth  rate  of  each  category.  The 
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total  growth  rate  is  then  calculated  and  stored  in  FHEFF. 

Output 

The  environemntal  parameters  listed 
under  the  Tables  and  Items  section  are  output  from  HEAT. 

Interfaces 


Subroutine  HEAT  interfaces  with  two 
main  subroutines,  INITIAL  and  BUDGET.  Subroutine  INITIAL 
is  used  to  retrieve  data  from  the  CRANDIC  data  base  and 
BUDGET  is  used  to  calculate  the  thin  and  thick  ice 
growth.  The  parameter  list  for  INITIAL  is  defined  as 
follows ; 


i) 

parameter 

1  -  number,  indicating 
which  catalog  name 

is  to  be  used  in 

INITIAL; 

ii) 

parameter 

2  -  array  used  to  store 

the  environemntal 

data  returned  from 

INITIAL; 

iii) 

parameter 

3  -  i  grid  points  of  the 
grid; 

iv) 

parameter 

4  -  j  grid  points  of  the 
grid; 

v) 

parameter 

5  -  DTG  required; 

vi) 

parameter 

6  -  number  of  points  in 
the  grid; 

vii) 

parameter 

7  -  tau  value. 
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as  follows; 


The  parameter  list  of  BUDGET  is  defined 


i) 

parameter 

1  - 

ii) 

parameter 

2  - 

iii) 

parameter 

3  - 

iv) 

parameter 

4  - 

iv) 

parameter 

5  - 

v) 

parameter 

6  - 

vi) 

parameter 

7  - 

vii) 

parameter 

8  - 

vii  i ) 

parameter 

9  - 

ix) 

parameter 

lO- 

x) 

parameter 

ll- 

Tables 

and  Items 

ice  thickness; 
growth  rate  returned 
to  HEAT; 

flag  used  to  de¬ 
termine  whether  thin 
or  thick  ice  growth 
is  calculated; 
grid  size; 
grid  size; 
wind  data; 
ice  temperature; 
mixed  layer  tempera¬ 
ture  ; 

air  temperature; 
surface  moisture; 
incoming  longwave. 


defined  in  HEAT; 


The  following  new  major  variables  are 


i) 

ii) 

iii) 

iv) 

v) 

vi) 

vii) 


vii  i) 


TICE  -  ice  temperature; 

TMI X  -  mixed  layer  temperature; 

TAIR  -  air  temperature; 

QA  -  surface  moisture; 

FLO  -  incoming  longwave 
radiation; 

PS  -  surface  pressure  and  sen¬ 
sible  heat  flux; 

CS  -  surface  vapor  pressure  and 
sensible  plus  evaporative 
heat  flux; 

UG  -  wind  values. 
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The  following  common  block  is  defined: 

/RAD/  -  contains  incoming  shortwave 
radiation. 


2.4.3.10  Subroutine  BUDGET 

Subroutine  BUDGET  is  used  to  calculate 
the  thin  and  thick  ice  growth  rates  by  using  a  simple  heat 
budget. 

Input 

Subroutine  BUDGET  receives  all  major 
input  varibales  from  HEAT.  Most  of  these  are  passed 
through  the  formal  parameter  list.  One  variable,  FSH ,  is 
passed  in  common  block  /RAD/. 

Processing 

The  initial  code  of  BUDGET  is  dedicated 
to  defining  all  necessary  constants  used  in  the  heat 
budget  calculations. 

A  branch  is  made,  depending  upon  the 
value  of  KCPEN.  if  KOPEN  is  less  than  zero,  processing 
continues  to  compute  the  thin  ice  growth  rate.  If  KOPEN 
is  greater  than  zero  processing  jumps  to  compute  the  thick 
ice  growth  rate. 

Continuing  with  the  the  ice  growth,  the 
main  heat  budget  equation  components,  and  growth  rate  de¬ 
rived  within  loop  101.  Subroutine  BUDGET  ends  processing 
here  for  this  branch. 


When  the  thick  ice  growth  rate  is  com¬ 
puted,  there  are  two  main  components  of  the  heat  budget 
calculation.  The  ice  temperature.  TICE,  is  solved  for 
iteratively  (Newton-Raplson  technique)  and  used  in  the 
head  budget  equation.  When  the  ice  temperature  values  are 
relatively  stable  between  iterations,  processing  finishes 
returning  the  thick  ice  growth  rate. 

Output 

The  output  of  BUDGET  is  contained  in 
the  following  variables; 

i)  FHEFF  -  thick  ice  growth  rate; 

ii)  FO  -  thin  ice  growth  rate. 

Interfaces 


routines. 


Subroutine  BUDGET  calls  no  other  sub¬ 


Tables  and  Items 


The  variable,  ALB,  defining  the  surface 
albedo,  is  the  only  major  variable  not  previously  defined 
and  used  within  BUDGET. 


2.4.3.11  Subroutine  GROWTH 

Subroutine  GROWTH  is  used  to  calculate 
the  change  in  ice  thickness  amd  concentration  due  to 
growth/decay. 


48 


All  required  input  variables  to  GROWTH 
are  passed  from  ICEMDL  through  the  formal  parameter  list. 
No  variables,  not  previously  defined,  are  input  to 
GROWTH. 

Processing 

The  processing  of  GROWTH  is  contained 
in  one  major  loop.  The  amount  of  ice  grown  and  melted 
during  the  time  step.  The  changes,  reflected  by  the 
growth/decay  rates,  are  added  to  the  ice  thickness,  HEFF, 
and  ice  concentration,  AREA,  variables.  The  ice  concen¬ 
tration  value  is  then  checked  to  contain  it  within  the 
limits  specified.  Ice  concentrations  are  not  allowed  to 
be  larger  than  1.0. 


Output 

The  variables  HEFF,  AREA  and  GAREA  con 
tain  output  variables  of  GROWTH. 

Interfaces 


GROWTH. 


No  other  subroutines _ -are  called  by 


Tables  and  Items 


GROWTH. 


No  new  variables  are  introduced  by 
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2.4.3.12 


Subroutine  ADJUST 


Subroutine  ADJUST  is  used  to  set  up  the 
thickness  and  concentration  at  the  outflow  cells. 

Input 

All  input  variables,  to  ADJUST  are 
passed  from  ICEMDL  through  the  parameter  list.  No 
variables  not  previously  defined  are  used. 

Processing 

Subroutine  ADJUST  is  called  at  the  end 
of  each  time  step.  ADJUST  uses  routine  MEAN  to  calculate 
the  ice  in  the  open  cells  by  taking  an  average  of  all  grid 
cells  adjacent  to  the  open  boundary.  The  ice  thickness 
and  concentration  arrays  are  both  modified  in  this  manner. 

Outputs 

The  arrays  HEFF  and  AREA  (ice  thickness 
and  ice  concentration)  are  respectively  modified. 

Interfaces 


Subroutine  ADJUST  calls  subroutine  MEAN 
to  calculate  the  ice  in  open  grid  cells. 

Tables  and  Items 


No  major  new  variables  are  defined. 
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2.4.3.13  Subroutine  MEAN 

Subroutine  MEAN  is  used  to  calculate 
the  amount  of  ice  in  open  cells. 

Input 

All  data  input  to  MEAN  is  passed  by 
ADJUST  through  the  parameter  list.  No  new  variables  are 
introduced  as  input. 


Processing 

Subroutine  MEAN  calculates  the  amount 
of  ice  in  a  grid  cell  as  the  mean  of  ice  in  the  adjacent 
cells.  Array  CUT,  which  is  the  boundary  mask  specifying 
the  outflow  cells,  is  used  to  control  the  calculations  to 
output  the  mean  ice  held  within  the  outflow  cells. 

Output 

The  variable,  HMEAN,  is  output  to 
ADJUST,  containing  the  mean  ice  thickness  and  concentra¬ 
tion  on  the  open  cells. 

Interfaces 


No  other  subroutines  are  called  by  MEAN. 
Tables  and  Items 

No  new  major  variables  are  introduced 

by  MEAN. 
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2.4.4 


Output  Module 


2. 4. 4.1 


Subroutine  DIVEKG 


Subroutine  DIVERG  is  used  to  output  the 
divergence  values  to  the  CRANDIO  output  file. 

Input 

All  required  input  variables  are 
supplied  to  DIVERG  from  ICEMDL  through  the  formal  para¬ 

meter  list.  No  new  variables  are  input. 

Processing 

Subroutine  DIVERG  operates  like  all 
output  module  subroutines.  Major  processing  functions  are 
performed  for  the  sole  purpose  of  setting  up  the  data  and 
appropriate  record  names  for  the  CRANDIG  data  base. 

The  data  is  placed  into  the  common 

block  array  DIVRG.  The  record  label  is  set  up  using  a 

data  statement  which  specifies  the  catalog  name,  tau  value 
and  record  length. 

If  an  error  occurs  on  the  CWRITER 

function  an  error  message  is  written,  stating  this  fact 
and  the  program  will  terminate. 

Output 

Subroutine  DIVERG  places  the  divergence 
values  on  the  CRANDIO  output  file. 
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Interfaces 


Subroutine  DIVERG  interfaces  with  the 
standard  FNOC  CRANDIC  software. 

DIVERG  is  called  on  every  fourth  time 

step  by  ICEMDL. 


Tables  and  Items 


The  following  common  block  is  defined; 


/DV/  -  MDV 

DIVRG 

FILLV 


contains  CRANDIO  ID  block; 
contains  data  values; 
fills  the  small  page  (requirement  of 
CRANDIC) . 


2. 4. 4. 2  Subroutine  UVPLOT 

Subroutine  UVPLOT  is  used  to  output  the 
ice  drift  forecasts. 


Input 

All  input  to  UVPLCT  is  provided  by 
ICEMDL,  through  the  parameter  list.  No  new  variables  are 
introduced . 

Processing 

The  first  main  function  of  UVPLOT  is  to 
convert  the  U,  V  ice  drift  components  to  ice  drift  direc¬ 
tion  and  speed.  This  is  accomplished  through  the  FNOC 
routine  UVDDFF .  Ice  drift  speed  is  converted  from  m/sec 
to  knots. 
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The  final  function  is  setting  up  of  the 
CRANDIC  record  lables  for  CRANDIC.  The  record  labels  are 
set  up  using  a  data  statement  containing  the  catalog 
names,  tau  values  and  record  lengths. 

If  an  error  occurs  on  the  CWRITER 
function,  a  message  is  written  and  the  program  terminates. 

Output 

The  ice  drift  direction  and  speed  is 
output  to  the  final  CRANDIO  file. 

Interfaces 


UVPLCT  interfaces  with  the  FNGC  CRANDIC 
software.  UVPLCT  is  called  on  every  fourth  time  step  by 
ICEMDL. 


Tables  and  Items 


by  UVPLOT. 


No  new  major  variables  are  introduced 


2.4.4. 3 


Subroutine  BUOYD 


Subroutine  BUOYD  is  used  to  calculate 
the  drift  of  simulated  buoys  placed  in  the  model. 

Input 


All  input  required  for  BUOYD  is 
supplied  by  ICEMDL  through  the  parameter  list.  No  new 
variables  are  introduced  by  BUOYD. 
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The  buoy  drift  positions  are  input 
through  common  block  /BUOY/ . 

Processing 

The  initial  processing  of  BUCYD  checks 
the  position  of  each  buoy.  This  is  done  to  determine  if  a 
buoy  has  moved  out  of  the  grid.  Once  a  buoy  reaches  a 
boundary  it  is  no  longer  tracked.  The  routine  branches 
around  the  remainder  of  the  processing  and  goes  on  to  a- 
nother  buoy  when  one  has  moved  to  the  boundary. 

Subroutine  INTRP,  a  FNOC  utility 
routine,  is  used  to  interpolate  the  U,  v  ice  drift  com¬ 
ponents  to  the  buoy  position  recorded  during  the  previous 
time  step.  The  u,  v  values  are  then  used  to  calculate 
how  far  the  buoy  shall  have  moved  during  the  current  time 
step.  This  distance  is  then  determined  in  terms  of  grid 
units  on  both  the  model  grid  and  the  FNOC  hemispheric  grid. 

The  final  position  of  each  buoy,  in 
terms  of  both  grid,  is  printed  in  a  table. 

Output 

The  positions  of  the  buoys,  calculated 
in  common  block  /BUOY/  are  output  from  BUOYD  a  printed 
table,  outlining  the  position  of  each  buoy  is  provided. 

Interface 

BUOYD  utilizes  the  FNOC  utility  routine 
INTRP  to  provide  u,  v  components  of  ice  drift  at  buoy  po¬ 
sitions. 


Tables  and  Interfaces 


BUOYD. 


No  new  variables  are  introduced  within 


2. 4. 4. 4 


Subroutine  HAPLOT 


Subroutine  HAPLOT  is  used  to  output  the 
ice  thickness  and  concentration  values  to  the  CRANDIC  file. 

Input 

All  required  input  to  HAPLOT  is  pro¬ 
vided  by  ICEMDL  through  the  formal  parameter  list.  No  new 
input  variables  are  defined. 

Processing 

Processing  within  HAPLOT  performs  two 
functions.  The  first  is  a  unit  conversion  while  the 
second  creates  the  necessary  data  for  CRANDIC.  Ice  thick¬ 
ness  is  output  first.  The  thickness  values  are  converted 
from  meters  to  cm.  The  CRANDIC  record  labels  and  data  are 
set  up  into  their  respective  positions  and  CVvRITER  is  used 
to  write  the  data. 


The  ice  concentration  values  are 
handled  in  the  same  manner.  However  no  unit  conversion  is 
performed . 

If  an  error  results  in  the  processing 
of  any  CWRITER  an  appropriate  message  is  written  and  the 
program  shall  terminate  processing. 
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Output 


The  ice  thickness  and  compactness 
values  are  output  to  the  CRANDIO  file. 

Interfaces 


software . 


HAPLOT  interfaces  with  the  CRANDIO 


Tables  and  Items 


No  new  major  variables  are  introduced. 


2 . 4 . 4 . 5  Subroutine  GROWDEC 


Subroutine  GROWDEC  is  used  to  output 
the  open  water  growth  and  ice  growth  forecasts. 

Input 

All  required  input  is  passed  to  GROWDEC 

by  ICEMDL. 


Processing 

The  processing  of  GRCWDEC  proceeds 
exactly  as  other  output  module  subroutines.  The  CRANDIO 
labels  are  defined  and  CWRITER  is  used  to  output  the  data. 

Output 

Forecasts  of  open  water  and  total  ice 
growth  are  written  to  the  CRANDIO  file. 
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Interfaces 


software. 


GROtoDEC  interfaces  with  the  CRANDIO 


Tables  and  Items 


Nom  new  variables  are  defined. 


2. 4. 4. 6  Subroutine  PRNT 

Subroutine  PRNT  is  a  small  subroutine 
which  is  used  to  print  model  arrays.  The  processing  of 
PRNT  depends  entirely  upon  the  input  data  specifications. 

Input 

The  formal  parameters  are  defined  as 

follows ; 


i) 

ARRAY 

array  to  be  printed; 

ii) 

I,J,K 

dimensions  of  ARRAY; 

iii) 

MI,MZ 

columns  of  ARRAY  which 

are  printed; 

iv) 

N 

-  number  of  rows  of 

ARRAY  to  be  printed. 

Output 

Printed  output  of  a  data  array  is  pro¬ 
vided. 
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Section  3.0 


Environment 


3.1  Equipment  Environment 

FNCC  operates  a  multiprogramming/- 
multi--  mainframe  computer  system  consisting  of  three  CDC 
6500 's,  with  131K  each  of  central  memory  (CM),  9  CDC  CYBER 
170/175  with  196K  CM  and  1  million  words  of  CDC  7030  ex¬ 
tended  core  storage  (ECS),  related  auxiliary  equipment  (7 
and  9  track  tape  units,  disk  storage,  etc.),  a  CYBER  203 
with  2  million  words  of  CM  and  its  front-end  processor 
CYBER  170/720,  and  a  VARIAN  plotter  and  its  plotting  soft¬ 
ware.  The  Sea  Ice  model  is  designed  to  run  on  the  CYBER 
203  computer,  and  its  output  is  plotted  on  the  VARIAN 
plotter . 


3.2  Support  Software 

FNOC  operates  under  the  NOS/BE  operat¬ 
ing  system.  This  system  contains  many  local  enhance¬ 
ments/modifications  to  facilitate  ease  of  operation.  Most 
of  the  enhancements  are  documented  in  either  the  FNOC  sub¬ 
routine/utility  file  or  the  FNWC  Computer  User  Guide 
Edition  2.  The  Sea  Ice  Model  was  converted  to  CDC  CYBER 
200  FORTRAN  Language,  version  1.5  and  utilizes  various 
routines  available  on  FNWCLIB. 


3.3 


Data  Base 


A  number  of  data  bases  are  maintained 
and  needed  by  the  sea  ice  model. 
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3.  3.1 


General  Characteristics 


Three  separate  input  data  bases  are 
required  for  the  sea  ice  model  to  properly  execute.  Two  of  these 
are  created  by  the  user.  The  third  data  base  consists  of  the 
FNOC  environmental  data.  Currently  this  data  is  not 

operationally  available  on  the  CY203,  therefore,  this  data  base 
is  also  set  up  by  the  user,  from  data  available  on  other 
machines. 


i)  TAPE7 

TAPE 7  is  the  mnemonic  for  the  main  input 
file.  This  file  is  set  up  by  a  user  and  remains  permanent  for 
the  desired  configuration  of  the  model  runs. 

The  file  is  utilized  by  the  input  module 
and  describes  the  execution  configuration  for  the  current  run 
(e.g.,  grid  size,  grid  location,  grid  configuration). 

ii)  TAPE8 

TAPE 8  is  a  more  dynamic  data  file  than 
TAPE7.  TAPE8  contains  the  DTG's  of  the  days  during  the  current 
run.  Therefore,  this  file  will  change  for  each  run. 

The  file,  TAPE8 ,  needs  to  be  defined  by 
the  user  of  the  sea  ice  model. 

iii)  FNOC  Data 

FNOC  data  is  kept  on  a  CRANDIO  file 
which,  at  the  time  of  this  writing,  must  be  set  up  by  the  user. 
These  data  will  change  for  a  specific  run.  The  data  is  obtained 
from  another  machine  which  has  access  to  the  FNOC  master  data 
base  (MASFNWC ) . 
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iv)  Output  Data-Base 


A  CRANDIC  output 

data 

base  is 

main- 

tained  by  the 

sea  ice  model.  This 

data 

base  is 

filled 

with  forecast 

variables  computed  by 

the 

sea  ice 

model. 

Variables  are  written  in  specified  time  steps,  defined  in 
routine  ICEMDL. 


3.3.2  Organization  and  Detailed  Description 

i)  TAPE7 

The  file  labled  TAPE7  is  a  binary  file 
accessed  by  a  formatted  READ.  The  following  information 
is  contained,  with  formats; 

i)  grid  size  specifications  -  215; 

ii)  FNOC  i  grid  points  for  MESH  - 
2F10.2; 

iii)  FNOC  j  grid  points  for  MESH  - 
2F10.2; 

iv)  number  of  rows/columns  in  model 
grid  -  15 

v)  uv  boundary  mask  -  27F2.0,  27  rows 

vi)  thickness  boundary  mask  -  28F2.0, 
28  rows; 

vii)  outflow  boundary  mask  -  28F2.0, 

28  rows; 

viii)  ocean  current  direction  -  25F3.0, 
25  rows. 

ii)  TAPE8 

The  file  labeled,  TAPE8,  is  a  binary 
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file  accessed  with  a  formatted  READ.  The  following  in¬ 
formation  is  contained,  with  formats; 

i)  number  of  time  steps  to  be  run  - 
I5; 

ii)  DTG ,  one  row  for  each  time  step 
3(A8,2X) . 

The  DTG  rows  are  defined  as  described 

under  1CEMDL. 


iii)  FNOC  Input  Data 


The  input  FNOC  data  base  is  a  CRANDIC 
data  file  on  the  CY203.  CRANDIO  is  the  operational  data 
file  format,  specified  by  FNOC,  on  the  CY203.  Detailed 
characteristics  of  the  CRANDIO  files  can  be  found  in  the 
CRANDIO  software  documentation  distributed  by  FNOC. 

The  data  on  the  CRANDIC  file  is  created 
by  transferring  ZRANDIC  data,  from  other  FNCC  machines  to 
the  CY203  with  proper  ZRANDIO  to  CRANDIO  conversion  soft¬ 
ware. 

The  following  records  are  contained  on 
the  CRANDIO  input  file; 


Record 


Contents 


A20 

A21 

A10 

A01 

A12 

All 

A18 

A16 


u-wind  component 

v  wind  component 

air  temperature 

surface  pressure 

surface  vapor  pressure 

short  wave  radiation  incoming 

total  heat  flux 

sensible  plus  evaporative  heat  flux. 
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iv)  CRANDIO  Output  File 

The  CRANDIO  output  file  is  created  by 
the  sea  ice  model.  CRANDIO  records  are  written  consisting 
of  forecast  variables  on  certain  time  steps.  The  follow¬ 
ing  records  are  written  every  four  time  steps; 


Record  Contents 


i) 

DIV 

ice 

divergence/convergence 

ii) 

FFF 

ice 

drift  speed 

iii) 

DDD 

ice 

drift  direction 

iv) 

THK 

ice 

thickness 

v) 

COM 

ice 

concentration 

vi) 

PRS 

ice 

pressure 

vii) 

GAR 

open  water  growth 

viii) 

HDF 

ice 

growth. 

Section  4.0 


Program  Maintenance  Procedures 


4 . 1  Conventions 

The  Sea  Ice  Model  system  adheres  to 
structured  design  and  programming  principles.  Flowcharting  and 
naming  conventions  adhere  to  the  standard  identified  below. 

a.  FNWC  User  Guide,  Edition  2, 

February,  1974. 

b.  CDC  Programming  Standards,  CDC-STD 
1.80.000,  December,  1971. 


4.2 


Verification  Procedures 


The  methods  of  verifying  the  sea  ice  mode.', 
output  through  display  of  the  output  file  on  plotting  display 
o  printed  output.  Plotting  of  ice  drift,  ice  growth/decay 
rate,  ice  thickness  is  a  very  efficient  method  used  to  check 
output  of  the  model. 


4.3 


Error  Conditions 


This  section  describes  the  error  conditions 
determined  by  the  Sea  Ice  Model. 


ICEMDL 

Message 

-  "OPEN 

error" 

Reason 

error 

in  opening  the 

CRANDIC  output  file 
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RELAX 


Message 

Reason 

UVPLOT 

Message 

Reason 

Result 

HAPLOT 

Message 

Reason 

Result 

DIVERG 

Message 

Reason 

Result 


"No  convergence  after  1300 
iterations" 

UICE  and  VICE  are  not  conver¬ 
gent  after  1300  iterations. 


"STATUS  is  (value)  ON  WRITE  CF 
filename" 

ISTAT  is  not  equal  to  0. 

Output  of  UVPLOT  is  not 
written  to  CRANDIO  file, 
program  stops. 


"STATUS  IS  (value)  ON  WRITE  OF 
filename" 

ISTAT  is  not  going  to  0. 

Output  of  HAPLCT  is  not 
written  to  CRANDIO  file,  pro¬ 
gram  stops. 


"CWRITE  STATUS  IS  (Value)  ON 
WRITE  OF  filename" 

ISTAT  is  not  equal  to  0,  the 
CRANDIO  output  file  is  not 
opened . 

Output  of  DIVERG  is  not 
written  to  CRANDIO  file,  pro¬ 
gram  stops. 


INITIAL 


Message  - 
Reason 

Result 


"CHECKNC  no  data  initial" 
No  required  data  field  for 
input 

Stop  the  program. 


4.4 


Special  Maintenance  Procedures 


There  are  no  special  maintenance  pro¬ 
cedures  for  the  Sea  Ice  Model  program. 


4 . 5  Special  Maintennance  Programs 

There  are  no  special  maintenance  pro¬ 
grams  for  the  Sea  Ice  Model. 


4 . 6  Listings 

Listings  of  the  Sea  Ice  Model  program 
and  subroutines  are  to  be  found  in  Appendix  C. 
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-:"****&whe 


Call 

GROWTH 


FELLIP 


> - - 


- 

/ _ 

Compute 

finite 

differences 

4 

i - ^ 

C  ena  ) 


Finite  differences 
for  RELAX  as  a 
function  of  current 
u.v. 


Q  FELLD1 


Compute 

finite 

differences 


Finite  differences 
for  main  interaction 
process  in  RELAX 


RELAX 


V _  V 

-  V 

Calculate 

Matrix 

Coefficients 


ur. 


Call 

FELLIP 

L _ 

Calculate 
variables 
of  previous 

U.V. 

L 

n 

Call 

FELLDI 

I 


Calculate 

new 

u.v. 


2 


DIVERG 


INITIAL 


Set  up 
CRANDIO 
record 
label 


"record 


Call 

INTRP 


HAPLOT 


y 


^  Buoy 
.off  gri 


1*1  If  Vi 


UNCLASSIFIED 

2o#2 


CONTROL  DATA  COUP  MONTEREY  CA  F/0  0/12 

PR09RAM  MAINTENANCE  MANUAL  POLAR  ICE  FORECAST  SUBSYSTEM  -  ARCTJ— CTC(U) 
OCT  81  PA  HARR#  T  C  PHAM#  J  P  WELSH  NOOO14-01-F-OO20 

NORDA-TN-122  NL 


END 


Appendix  B 


Grid  Structure 


The  model  grid  is  set  up  in  a  staggered  manner.  Momen¬ 
tum  variables  are  defined  at  the  grid  points  while  thermody¬ 
namic  variables  are  defined  for  grid  cells. 

Environmental  variables  accessed  from  the  FNOC  data  base 
are  computed  to  be  valid  at  the  grid  cell  locations.  Subvent- 
ive  AVG  is  used  to  compute  the  grid  cell  averages  for  the 
thermodynamic  variables. 

The  following  example  illustrates  the  definition  of 
grid  parameters  for  a  small  grid. 

Say,  we  define  NX, NY  to  be  7  (dimension  of  momentum  var¬ 
iable  grid).  Figure  Bl  illustrates  that  NX1,NY1  (thermody¬ 
namic  grid)  will  be  8  while  the  value  of  NUMBER  will  be  9. 
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Figure  Bl  Grid  Configuration 
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